Field-dependence of magnon decay in yttrium iron garnet thin films 
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We discuss threshold field-dependence of the decay rate of the uniform magnon mode in yttrium 
iron garnet (YIG) thin films. We demonstrate that decays must cease to exist in YIG films of 
thickness less than 1 (im, the lengthscale defined by the exchange length. We show that due to the 
symmetry of the three-magnon coupling the decay rate is linear in AH — (H c — H) in the vicinity 
of the threshold field H c instead of the step-like Fk=o ocO (AH) expected from the two-dimensional 
character of magnon excitations in such films. For thicker films, the decay rate should exhibit 
multiple steps due to thresholds for decays into a sequence of the two-dimensional magnon bands. 
For yet thicker films, such thresholds merge and crossover to the three-dimensional single-mode 
behavior: r k=0 oc |A#| 3/2 . 



PACS numbers: 75.10.Jm, 75.40.Gb, 75.30.Ds 

Introduction. — Extensive experimental and theoretical 
research on the ferromagnetic insulator, yttrium iron gar- 
net [Y3Fe2(Fe04)3 or YIG] was started more than half 
a century ago has benefited from this material's excep- 
tional purity, high Curie temperature, and relative sim- 
plicity of the low-energy magnon spectrum.^ Recent dis- 
covery of the Bose-Einstein condensation of the highly 
occupied magnon states created by microwave pumping 
in YIG thin films has attracted substantial interest P 
Recently, threshold effects due to the so-called three- 
magnon splitting have been reported in a quasi-two- 
dimensional (2D) thin films of YIG under microwave 
pumping and as a function of external fieldP In this, as 
well as in the other recent worksj^l control over the spin 
current enhancement in layered metal-ferromagnet struc- 
tures by the three-magnon processes is sought. This calls 
for a deeper theoretical insight into the decay dynamics 
of such structures. Fundamentally, given its outstand- 
ing properties,^ YIG may offer a fertile playground in 
the studies of threshold phenomena, because the decay 
conditions in it can be varied continuously by both film 
thickness and external magnetic field. 

Decays. — In this work, we discuss the decay rate of the 
uniform mode (k — magnon) in an insulating ferromag- 
netic thin film as a function of external magnetic field and 
film thickness. In particular, using magnon dispersion in 
the lowest-mode approximation, we outline the ranges 
of the field and film thickness that favor decays in YIG. 
Specifically, we show that kinematic conditions for three- 
magnon decays cannot be met for the YIG films of thick- 
ness d m i n fa 1 fj,m or less and for the fields iJ™ ax ps 600 Oe 
or higher. The upper limit on the external field is defined 
solely by the magnetization of a ferromagnet, in agree- 
ment with Ref. |21 Less obvious is the existence of the 
limit on the film thickness, which is fixed by another fun- 
damental characteristics of a ferromagnet: its exchange 
length. The physical reason for that limit is the decreas- 
ing role of the long-range dipolar interactions with the 
decrease of the film thickness. The presence of such a 
limit must be important for the control of relaxation and 
transfer of the spin current in layered structures, which 



rely on the three-magnon processes in YIG. 4 

We find that the threshold field-dependence of the de- 
cay rate for the k = magnon near the threshold field 
H c is r k=0 oc \AH\ ■ Q(AH), where AH = (H c - H), 
because the three-magnon interaction vanishes along the 
direction of the film magnetization, which is also precisely 
the k-direction where magnon band minima are located. 
This leads to a reduction of the phase space for decays 
and results in a vanishing decay rate near the threshold, 
contrary the naive expectation of the step-like increase of 
Fk=o oc@(AH), when the symmetry of the three-magnon 
interaction is ignored.^ We have supported our considera- 
tion by explicit calculations of both T = and T = 300K 
relaxation rate dependencies on the field for several rep- 
resentative YIG film thicknesses. 

The finite-size quantization in the film thickness direc- 
tion leads to the formation of multiple magnon bands. 6 
We argue that the field-dependence for the decay rate 
should exhibit multiple steps linear in \H Ci —H\, corre- 
sponding to the H Ci thresholds for decays into a sequence 
of magnon bands. For thick films, these multiple thresh- 
olds should merge in a continuum, and the decay rate will 
crossover to a three-dimensional (3D) single-band behav- 
ior: Fk=o oc |Aif| 3 / 2 . These results should be helpful 
in finding an optimal set of parameters for spin-current 
enhancement ™ 

Dispersion, density of states. — Despite its fairly com- 
plicated crystal structure, at low energies YIG can be 
described with great success as an effective large-spin 
Heisenberg ferromagnet on a cubic lattice with nearest- 
neighbor exchange, long-r ang e dipolar interactions, and 
negligible spin anisotropyP^ Thus, at long wavelength, 
magnon energy in YIG is determined by a competition 
between three couplings: exchange, dipolar, and Zeeman. 
For a ferromagnetic crystal of the film geometry and 
external field directed in-plane where it co-aligns with 
the magnetization direction, as is done most commonly 
in experiments, the lowest-mode approximation for the 
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FIG. 1. (Color online) (a) 2D plot of magnon dispersion in 
YIG in the lowest-mode approximation, Eq. for H — 1000 
Oe, d = 5 (im, and in k z > sector. Inset: axes conventions 
relative to the film and the field/magnetization direction, (b) 
2D magnon density of states (DoS) for several representative 
field values, E is in GHz. 
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FIG. 2. (Color online) d — H decay diagram for YIG, shaded 
area is where decays are allowed. Solid (dashed) boundary is 
the numerical (analytical) solution for the threshold bound- 
ary, see text. H max and dmin are shown. Inset: 2E ul i n /Eo vs 
H for several d's, upper threshold fields are indicated. 



magnon energy yields EEISI2I 

£ k = JJh + pk 2 + A k sin 2 k ) (h + pk 2 + A k ) , (1) 

where A k = / k A and A k = (1 — / k )A and h = pH, p — 
J So 2 , and A = 4irp,M are the energy scale parametriza- 
tions of the external field, exchange, and dipolar interac- 
tions, respectively. The form factor / k = (1 — e _ l k l d )/|k|c? 
is from the dipolar sums in the direction of the film thick- 
ness d 3> a, and # k is the angle between the ferromagnet's 
magnetization (directed in-plane, z axis) and magnon's 
in-plane 2D wave vector k, see Fig.[lja). We use p — gps 
where g = 2 is an effective g factor and pb is the Bohr 
magneton. In this work we adhere to the notations and 
units of Ref. [6j which has provided a detailed micro- 
scopic spin- wave theory of YIG in 1/S approximation, 
and we use experimental parameters for YIG, magnetiza- 
tion 4nM = 1750 G, exchange stiffness p/p = 5.17- 10" 13 
Oc m 2 , and lattice constant a = 12.376 A. In Fig. [TJa), 
the magnon dispersion from Eq. ([T]) for a representative 
field H — 1000 Oe and film thickness d = 5 pm is shown. 

It should be noted that the dipolar interactions are re- 
sponsible for the nontrivial structure of magnon band 
with the minima at finite wave vectors, which allow 
the decay conditions to occur. Although approximate, 
Eq. ([!]) provides a close quantitative description of the 
lowest 2D magnon energy band,^ quantized due to fi- 
nite film thickness d in the x direction [Fig. [TJa) , inset] . 
At small k, dipolar interactions dominate over the ex- 
change and result in a steep decrease from the energy of 
the uniform mode, Eq = \Jhih + A), for k's along the 
magnetization direction. At larger k, exchange energy 
dominates, giving s» h + pk 2 : while at intermediate 
k, competition between exchange and dipolar terms re- 
sults in peculiarly shaped minima at ±k m = (0, ±fe m ), 
see Fig.ilja). 



The 2D density of magnon states from Eq. ([lj is shown 
in Fig. |ljb) for several representative field values. Pre- 
dictably, 2D DoS exhibits a step-like increase at the band 
minimum, which is followed by an unusual, almost linear 
increase, reminiscent of the similar behavior for the rela- 
tivistic dispersion. The latter is due to a nonparabolicity 
of the long- wavelength magnon dispersion, see Fig. [lja) . 

"Decay diagram". — For the decays to take place the 
kinematic conditions must be fulfilled. For the two- 
magnon decay (thrcc-magnon splitting) of the uniform 
mode, the condition to be satisfied is simply Eq = 2Eq. 
With the microscopic parameters, such as exchange stiff- 
ness and magnetization, fixed, some other parameters 
can be varied to allow or to forbid decays altogether. 
As is clear from Eq. ([T]), the external field increases the 
energies of both the uniform mode and the minimum, 
making decays kinematically impossible at some higher 
field valueP Another parameter is the film thickness d, 
which enters Eq. ([!]) through the form factor / k . While 
the manner in which d influences decays is not a priori 
clear, both trends, vs field and vs thickness, can be ex- 
amined numerically. Such an examination is exemplified 
in Fig. [2j which shows 2E m i n /Eo vs field for several film 
thicknesses. Clearly, when the plotted quantity is < 1, 
decays are allowed, and the crossing of 1 corresponds to 
a threshold field for decays. 

Our Fig.[2]gives the complete d—H "decay diagram" for 
YIG with the shaded area showing the parameter space 
where decays are allowed. Two key results are clear in 
Fig. [2} (i) the upper threshold field does not exceed some 
ffnmx even f or i ar g e values of d, and (ii) there exists a 
lower limit on the film thickness d mm , below which the 
decays of the uniform mode may not occur at all. The 
solid boundary is the numerical solution of Eq. ([T]) for the 
energy minimum and the decay conditions. The dashed 
line is an approximate analytical solution, which turns 
out to be very precise. The latter also gives us a deeper 
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insight into the nature of 7J™ ax and d m - ln discussed next. 

Large-\k m \d approximation — At large enough d the 
wavevector of the magnon energy minimum satisfies 
|k m |c!^> 1. Then the formfactor /k m ~ l/|k m |d is small, 
reflecting the reduced role of dipolar interactions in the 
energy of the k m -magnon. One can then show^ th a t. 
both the exchange and the dipolar energies for the km- 
magnon scale as oc d~ 2 ' 3 and thus are -C h for any reason- 
able field. This implies that the energy of the magnon 
band minimum is E mm w h, expected for the uniform 
mode in the absence of dipolar interactions. Then, the 
decay threshold equation, 2E mnl — Eq, trivially gives 
h c = A/3, relating saturated value of the decay thresh- 
old field to the material's magnetization: i?™ ax = 47rM/3 
(= 583 Oe for YIG in Fig. gb)), see also Ref. H 

The physical question is: what parameter of the ferro- 
magnet defines d m ; n ? One can extend the large- |k m |d 
approach and fincP^I that the energies of the dipolar 
and exchange interactions must be related by |k m |d = 

1/3 

(Ad 2 /Ap) . With this, the decay threshold condition 
2-E m in = Eq leads to an algebraic equation in H c vs d, 
which can be resolved in a compact form. 1 *-' It is plot- 
ted as a dashed line in Fig. [2jb), which coincides almost 
exactly with the decay boundary obtained from Eq. 
numerically. From the same equation we find the mini- 
mal thickness to be d m in~C-\/p/A, explicitly related to 
the exchange length of the ferromagnet, £ ex = A /p/A, al- 
beit with a large numerical coefficient C « 62. 04pSJ Using 
parameters for YIG, the exchange length is £ ex « 13. 9o 
and the minimum thickness <i m i n = 1.067 fim, remarkably 
close to the numerical result <i m i n = 1.017 /mi. 

The physical reason for the very existence of d m i n is the 
decreasing role of long-range dipolar interactions in the 
magnon's energy with the decrease of film thickness, as 
the relation of rf m i n to exchange length implies. The fact 
that the decay boundary in d exceeds the exchange length 
by a large numerical factor is due to the rather stringent 
requirements of the decay conditions. We emphasize that 
the provided d—H diagram should apply equally to the 
other thin-film ferromagnets. 

Decay rate. — Transitions that involve changing the 
number of magnons, such as decays, recombination, or 
coalescence, originate from the dipolar interactions that 
couple longitudinal and transvers e sp in components and 
therefore do not conserve magnetic^^ as well as mechani- 
cal angular momentumP Microscopically, dipolar interac- 
tions result in anharmonic couplings of magnons, which, 
for the decay processes of the k = uniform mode into 
two magnons at q and — q, can be written as 

^ 3) -^E<U(<«Vo+H.c). (2) 
q 

The three-magnon coupling in Eq. ^ has an angular de- 
pendence: V ( S,-q = Vosin26» q , with V = A/v^! 11111111 
This angular dependence is essential since it reflects the 
symmetry of dipolar coupling of the transverse, S x (S v ), 
and longitudinal, S z , spin components: V xz oc xz/r 3 
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FIG. 3. (Color online) The T = decay rate r k = vs H for 
d = 2, 5, and 10 /im. Dotted line is |Tk=o for d = 10 /jm 
with the angular-dependence of the three-magnon coupling 
omitted, V^'.q => V . 

(V yz oc yz/r 3 ). In particular, it is natural for this cou- 
pling to vanish for the spin-wave propagating with the 
momentum q along the direction of magnetization M [z 
axis, Fig. nJa)]pJ We would like to point out that it is 
also precisely the direction along which the minima of 
the magnon band are located.^ Therefore, the ampli- 
tude of the decay of k = magnon into two magnons at 
the band minima q TO and — q m is zero. This will lead to 
a rather spectacular violation of the naive expectation: 
while kinematic conditions for the decay of k = magnon 
into ±q m are just met at H c , the corresponding decay 
amplitude is vanishing. Thus, the decay rate must in- 
crease gradually from the threshold, not in a jump-like 
fashion as in the DoS, Fig. [jib). 

At T = only spontaneous magnon decays are 
allowedP The three -magnon recombination processes 
have to obey the same kinematic constraints as the de- 
cay, having therefore the same threshold conditions. The 
three-magnon coalescence processes involving the k = 
mode correspond to the "vertical" transitions, which are 
forbidden either kinematically as in the single mode case 
or by the quantum number of the interband transition 
in the multi-band situation. The four-magnon scattering 
amplitude from the exchange interaction vanishes iden- 
tically for the uniform mode, while the remaining four- 
magnon interactions from the dipole-dipole interaction 
together with impurity scattering should be providing 
a background with weak H and T dependence, distinct 
from the threshold behavior discussed here. 

With this in mind, using the kinetic approach}^ which 
takes into account the balance between decay and recom- 
bination processes in the relaxation time approximation, 
the decay rate^ of the uniform mode is 

r k=0 = tt J2 I ^oS-q 1 2 [ 2 "q + 1] 6 ( E o ~ 2£ q ) (3) 
q 

where n q = [e hEci ^ kbT — l]" 1 is the Bose occupation fac- 
tor. It is clear from Eq. (J3j) that while the magnitude of 
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FIG. 4. (Color online) Same as in Fig. [3] but for T = 300 K. 

the decay rate at finite T can be substantially modified 
from the T — result by the Bose-occupation factors, the 
qualitative threshold behavior must remain the same. 

One can investigate the threshold behavior of Eq. |3]) 
analytically and obtain Fk = ooc \Eq — 2£ l m j n | for H — > if c ™ 
Given the proportionality between (Eq — 2£ , m j n ) and 
AH =(H C — H) demonstrated in Fig. [2j this yields linear 
dependence of the decay rate on the field relative to the 
threshold: r k=0 oc |AF| • <d(AH), see Figs. [3] and [I] Once 
again, this is the consequence of an effective suppression 
of the phase space for decays due to the discussed angular 
dependence of the three-magnon coupling. 

In Figs. [3] and [4] we show T^ =0 for three different film 
thicknesses and for T = and T = 300 K, respectively. 
While the overall scale in the two figures is in a com- 
pletely different frequency range™ the shapes of T vs H 
are qualitatively very similar, especially concerning their 
threshold behavior vs field, in agreement with the above 
analysis. The relative difference between the curves for 
different d's reflects the smaller phase space for decays in 
thinner films. In the same Figs. [3] and [4] we demonstrate a 
dramatic contrast with the results of Eq. ^ if the angu- 
lar dependence of the three-magnon coupling in Eq. (pi is 

neglected, Vq ; =r-Vb. The latter results exhibit jumps 
at H^s and a linear increase after that, similar to the 
2D magnon DoS in Fig. [ljb). One should also note that 
the overall decay rate is also markedly overestimated if 
the angular dependence of the three-magnon interaction 



is ignored. 

In thicker films, the decay rate will be further modified 
by the multiple magnon bands that occur due to finite- 
size quantization in the film thickness direction.^ It will 
exhibit a sequence of steps linear in \H Ci —H\, where H Ci 
is the threshold field for decay into an ith band, increas- 
ing the decay rate every time the corresponding kine- 
matic conditions are met. Strictly speaking, the angle 

(3) 

9q in V . q _ q is between the 3D q vector and magnetiza- 
tion vector M. In the quasi-2D geometry with the levels 
quantized in the x directions, q x has discrete values and 
the minimal value of 0°i m 9™" /111 is zero only for 
the lowest magnon band. Thus, the thresholds in H Ci 's 
will have some small step-like behavior. However, this 
effect must be negligibly smaller compared to the steps 
111 Figs. |] and when the angular dependence in is 
ignored altogether. Extending our analysis to the limit 
of thicker films where multiple bands merge into a single 
3D band, using Eq. ^ we obtain^ that the decay rate 
near the threshold field will crossover to rk = o°c \AH\ 3 / 2 . 

Conclusions. — In this work, we have discussed the 
field-dependence of the decay rate of the uniform magnon 
mode in YIG thin films and the effects of film thickness 
in it. As a result of our analysis, we have established that 
the two key characteristics of a ferromagnet, its magneti- 
zation and exchange length, define the extent of the d—H 
parameter range that favors decays in thin film geome- 
try. Our calculations of the decay rate should provide an 
important guidance for the experimentalists in designing 
the optimal conditions for the control of spin current and 
its relaxation in thin films. A particularly intriguing sug- 
gestion to pursue is to study the properties of a thin film 
with varying thickness, which may permit decays and, as 
a consequence, the spin current enhancement in one part 
of the film and forbid it in the other. 
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I. SUPPLEMENTAL MATERIAL 

A. large- |k m |d approximation 

Assuming that for large enough d the magnon en- 
ergy minimum satisfies |k m |d ^> 1, makes the formfactor 
/ k = (1 - e-^ d )/\k m \d l/|k m |d small and allows 
to simply _E k in Eq. ([!]) to an algebraic expression in 
reduced energy scales of the dipolar and exchange inter- 
actions, A = A/|k m |d and p = p\\i m \ 2 : 



h + p)[h + p + A , 



(1) 



where we also used the fact that the minimum is located 
along the z-axis, so 6y m = 0. The minimum condition 
dE k /dk = yields: 



A u A 
n + p + — = n ■ — 

4 Ap 



(2) 



which, in the limit of both reduced exchange and dipolar 
energy scales being small compared to the field, p, A <C h, 
yields a straightforward relation between the two: 



A = 4p. 



(3) 



This relation gives the field-independent answer for 
|k TO |d: 



Ik™ Id = 



/Ad 2 
\4p 



2\ 1/3 



(4) 



which justifies the above approximation p, A <C h for 
large enough d because: 



A = Ap = 



4pA 2 
~d? 



2\ 1/3 



(5) 



Having the dipolar and exchange scales related by 
Eq. ([3]), we can substitute it into Eq. ([I]) for .Emm and 
use it to solve the decay threshold condition 2E m j n = E 
with Eq = \]h{Jn + A). This leads to a quadratic equa- 
tion in the threshold field h c (= pH c ) vs A (which, in 
turn, is a function of d by Eq. (|5|), that can be resolved 
in a compact form: 



G 




12 



(6) 



with two solutions being the upper and the lower thresh- 
old fields. They give the (approximate) threshold bound- 
ary in Fig. [2^b), shown by the dashed line. From the 
same solution one can find the "maximal" A m that cor- 
responds to the lower boundary on d = d mm for which 
the solution of the threshold condition exists. This is also 
the point at which h Cl = h C2 in Eq. ([6]) and the content 
of the square-root is zero, solving for which gives 



A TO = A(d min ) = 



6- V15 



21 



A. 



Using Eq. ^ we finally obtain: 

_ , , 21 \ 3/2 
6-V15 



62.04 , 



(7) 



(8) 



Although the following result can be trivially obtained 
from the solution for h ci in Eq. ^ taking the d — >• 00 
limit (A — > 0), one can derive the large-d limit for H Cl — > 
i?™ ax via an observation that at d — > 00 energy minimum 
in Eq. |TJ is E min w h, as A, p — > 0. This reduces the 
decay condition to a simple one: 



2h c = yj h c (h c + A), 



(9) 



with two solutions: h C2 = and h Cl = pH™ &x = A/3, 
thus relating the saturated magnetization of the material 
to the upper limit of the threshold field for decays. 



B. threshold field-dependence of Pk=o 

Since we are interested in the threshold behavior, only 
spontaneous (T — 0) decay rate is considered. Near the 
threshold, decays are happening into the magnons in the 
vicinity of the magnon band minima ±q m ~ (0, ±q m ). 
Expanding in (q — q m ) near the minima gives: 



2E n 



(10) 



2 



^ J dq y J dq z sin 2 26» q 5 yAE 



Aq 2 



where Aq = (q— q TO ), AE = (Eq — 2E min ) and m is the 
effective magnon mass. Rewriting sin20 q as 2g y g z /q 2 
and shifting integration in q z by q m yields 



rk=o = 



2^o 2 



7r|q« 



ql dq y dq z 8[ AE 



m 



(11) 
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FIG. 5. (Color online) A sketch of Fk=o vs H for the case of 
multiple magnon bands. Dotted lines show linear in \H Ci —H\ 
contributions of each band, solid line with the shading is the 
total effect, dashed line is the 3D limit oc AH 3 * 12 . 

where q z — q z — q n , q 2 = q 2 + q 2 , and in approximating 
q 2 ~ q 2 n the smallness of q z 's and q y 's in comparison 



with |q m | was used. It finally leads to: 
V 2 m 2 

Tk=o = n ry \Eq — 2£ , m ; n . (12) 

|q m | 

Since \E — 2E m i n \ oc \H C — H\, this yields r k=0 oc AH. 
At finite T ^> E min , the decay rate near the threshold 

is 

_ 2T 2 V 2 m 2 , , 

1 k=0 — i 12~ I _ ^min| ■ \i-Oj 

^min 11™ I 

For the 3D case, when the 2D bands merge into a con- 
tinuum, the above consideration results in: 

v ™^S?! qldq ^ 5 ^ E -^)> (l4) 

where q 2 = q 2 + q\ now and possible mass anisotropics, 
while neglected, are not going to change the result qual- 
itatively: 

r 3D „ Vom 2 | ,3/2 

1 k=0 K 1 ry | ^0 — ^ ft min | ) (-LOJ 

and thus, rg2 cx AiJ 3 / 2 . 



